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O ■ Abstract 

D ! 

The GR@PPA event generator has been updated to version 2.7. This distribu- 
tion provides event generators for V (W or Z) + jets (< 4 jets), VV + jets (< 2 
jets) and QCD multi-jet (< 4 jets) production processes at pp and pp collisions, in 
addition to the four bottom quark productions implemented in our previous work 
(GR@PPA_4b). Also included are the top-pair and top-pair + jet production pro- 
cesses, where the correlation between the decay products are fully reproduced at 
the tree level. Namely, processes up to seven-body productions can be simulated, 
based on ordinary Feynman diagram calculations at the tree level. In this version, 
the GR@PPA framework and the process dependent matrix-element routines are 
separately provided. This makes it easier to add further new processes, and allows 
users to make a choice of processes to implement. This version also has several new 
features to handle complicated multi-body production processes. A systematic way 
to combine many subprocesses to a single base-subprocess has been introduced, and 
a new method has been adopted to calculate the color factors of complicated QCD 
processes. They speed up the calculation significantly. 
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1 Introduction 



The great success of the Standard Model in recent decades leaves no doubt that gauge 
theories are capable of describing the interactions between elementary particles. Precise 
perturbative calculations based on gauge theories will become a crucial portion as larger 
collision energies are available to probe higher energy phenomena. For instance, the Higgs 
boson(s) and new particles predicted by theories beyond the Standard Model are the most 
important subjects to study in the current and future hadron collider experiments. They 
frequently decay to final states including multiple hadron jets, which predominant QCD 
interactions can easily imitate. A precise evaluation of these backgrounds is, therefore, 
necessary in order to accomplish a reliable measurement. However, once we try to do that 
using perturbative gauge theories, we immediately encounter a huge number of processes 
to evaluate. It is often too large to calculate by hand. 

We have been carrying on automatic computations of Feynman diagrams by using the 
GRACE system . GRACE is capable of calculations at the one- loop level [2 as well as 
at the tree level in the standard electroweak theory, and at the tree level in the minimal 
super symmetric extension of the Standard Model (MSSM) It can also be applied to 
QCD interactions. However, since the development has been mostly aimed at applications 
to lepton collisions, it is not directly applicable to hadron-collision interactions, mainly due 
to the necessity to include the parton distribution function (PDF). Besides, processes in 
hadron collisions usually consist of lots of subprocesses. They are desired to be combined 
into a single process. 

In order to implement those features specific to hadron collisions, we have developed an 
extended framework, called GR@PPA (GRace At PP/Anti-p). Early developments can 
be seen in our previous reports jUE]. The primary function of GR@PPA is to determine 
the initial and final state partons, i.e., their flavors and momenta of the incoming partons 
by referring to a PDF, and those of the final state partons for jets or decay products. 
Based on the GRACE output codes, GR@PPA calculates the cross section and generates 
unweighted parton-level events using BASES/SPRING [6 j. The GR@PPA framework also 
includes an interface to the LHA common data format jjj- The final steps of the event 
generation, i.e., the initial- and final-state radiation, hadronization, decay and so forth, 
are done by passing the unweighted partonic events to PYTHIA [8 or HERWIG jU]. 

Although the GR@PPA framework is not process-specific, we provide it as an event 
generator package including matrix elements for some selected processes. At the moment 
(version 2.7), the included processes are the weak boson(s) plus N jets and tt plus N jets 
productions as well as the pure QCD N jets productions, where N < A. These processes 
are the most important backgrounds in the Higgs boson and SUSY particle studies, and 
also useful to understand multi-body particle dynamics in precise measurements. Our 
previous work, the four bottom quark production (GR@PPA_4b [TU] ) , is also included. 

The reasons why we provide particular processes apart from the benefit of the auto- 
matic calculation by GRACE are the following; First, kinematical singularities in each 
process must be cared with a proper treatment. One can get a high efficiency to generate 
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unweighted events since the kinematics are already well optimized. This is immediately 
addressed to the program running speed. It is critical for large scale MC productions. 
Second, it will be easier to adopt higher order calculations. Once the customized matrix 
element for an NLO process is prepared [TT], users will be able to simply use it without 
any detailed care. Third, some different extensions are possible only by the modification 
of the framework. For example, GR@PPA generators can be extended to work under the 
C++ environment just by rewriting the framework in C++. The parton shower algo- 
rithm, e.g. the NLL parton shower will also be possible to implement by modifying 
the framework. 

We describe a symbolic treatment of the parton flavor in the diagram calculation in 
the next section. The feature of GR@PPA and its running are given in section 3. Some 
benchmark cross sections and program performances are presented in Sections 4 and 5, 
respectively. Finally, a summary is given in Section 6. We provide a list of supported 
processes and a summary of benchmark test conditions in Appendices. 

2 New features in the extension of GRACE to pp/pp 
collisions 

In hadron-hadron collisions, a certain process of interest usually consists of several in- 
coherent subprocesses. The total cross section is thus expressed as a summation of the 
subprocesses as 




where f?{x a , Q 2 ) is the PDF of the hadron a (p or p), which gives the probability to find 
the parton % with the energy fraction x a at the probing virtuality of Q 2 . The differential 
cross section d&ij-^F(s) j 'd&p describes the parton-level hard interaction producing the 
final-state F from the collision of partons, i and j, where s is the square of the total 
initial 4- momentum. The sum is taken over all relevant combinations of i, j and F. 

The original GRACE system assumes that both the initial and final states are well- 
defined. Hence, it can be applied to evaluating daij^F^s) / d^F and its integration over 
the final-state phase space only. An adequate extension is necessary to take into 
account the variation of the initial and final states. We have made two sorts of develop- 
ment in GR@PPA, — applying PDF in the phase space integration and sharing several 
subprocesses as a single base-subprocess. The former is described in our previous paper 
[TUj . Here, we focus on the later. As already mentioned, a "process" of interest is usu- 
ally composed of several incoherent subprocesses in hadron interactions. In many cases, 
the difference between the subprocesses is only in the quark combination in the initial 
and/or final states. The matrix element of these subprocesses is frequently identical, or 
the difference is only in a few coupling parameters and/or masses. In such cases, it is 
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convenient to add one more integration/differentiation variable to replace the summation 
in Eq. ((TJ with an integration. As a result, these subprocesses can share an identical 
" GRACE output code" and can be treated as a single base-subprocess. This modification 
drastically saves the program running time and simplifies the program coding. 

The number of combinations to take N out of M flavors, allowing an overlap, is given 
by mHn{= nWM^w )- I n case a ^ P ar t° n flavors up to the 6-quark are considered, 
M is equal to d, c, s, b, g, u, d, c, s, h). The configuration of the N jets final state has 
n-ffjv subprocesses. Clearly a smaller M decreases the number of subprocesses. Since the 
diagram structure is invariant against the quark (anti-quark) flavor exchange in QCD, 
the invariance is preserved by introducing generic up-type and down-type quarks even if 
the electroweak interaction is included. The base-subprocesses can be configured using 
these generic quarks. Then, the number of configurations is reduced to ^Hn <C nH^. 
The output code of the matrix element from the GRACE has been so modified that the 
masses and couplings can be treated as input variables. The number of combinations has 
been further reduced using symmetries in the Standard Model such as the invariance in 
the charge and parity conjugate. In Table^we list up the number of base-subprocesses for 
the N jets production process in pp(p) collisions, to be compared with the total number 
of subprocesses that they cover. The initial colliding partons in the base-subprocesses in 
Tabled are composed of q u q u (qdqd), q u qd, q u g{qdg), q u qd, q u q u {qdqd) and gg, where q u {qd) 
and g are the up(down)-type quark, and the gluon, respectively. 

Equation (J2J) can be rewritten in term of the base-subprocess using a weight factor to 
treat the initial and final state parton configuration as 

f f f - dalfZ C F d (s;rn,a) 

a — y / dx\ I dx 2 / d&p WijF — * , (2) 

f-tj J J d$ F 

where d<3fj^ ed is the differential cross section of the base-subprocess with input argu- 
ments of masses and couplings. The coefficient w^p is the weight factor for the chosen 
configuration. An appropriate subprocess is selected in event by event based on the flavor 
configuration decided by the weight factor. Since the QCD interactions are identical for 
any configuration in the subprocess, the weight factor is composed of PDF and the CKM 
(Cabibbo-Kobayashi-Maskawa) jT^j matrix as 

WijF = fl(x 1 ,Q 2 )ff(x 2 ,Q 2 )\VcKM\ 2K {Bi(X -> F') x r£ t } L , (3) 

where L is the number of X bosons (W or Z). The probability of the decay final state F' 
can be handled in terms of the branching ratio, if the interference with the other partons 
is ignored. The branching ratio and the total width can be given by the experimentally 
measured ones. The flavor configuration is determined by the |Vckm| 2K , where K is the 
number of W bosons in the process. 

We have to define a certain QCD color base representation in order to calculate the 
color factor of the interaction. In addition, PS programs usually require the information 
of color connection between partons in the event. The color bases should be so defined 
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that the required connection information can be easily derived. In our previous works 
[Hj we adopted such a definition that all color bases are composed of the SU(3) triplets. 
Namely, gluons are decomposed to a pair of quark and anti-quark. This definition includes 
color-singlet components which are not well interpreted for the hadronization. Besides, 
this leads to a large number of color bases for complicated processes consisting of many 
partons. 

In the present version of GR@PPA we have adopted a different definition. We use a 
chain of a quark, n gluons (n = 0, 1, 2, ...) and an anti-quark as fundamental components. 
The color bases are given by products of such chains, and "glueballs" having no quark at 
the ends. This definition fits to the requirement of the LHA event interface jjj- In the 
event generation, one of the possible color bases is chosen in proportion to their squared 
amplitudes in order to determine the color connection. The interference is ignored (i.e., 
a large N c approximation) as usual, since no reasonable way is known to include it. Of 
course, the interference is taken into account in the cross section evaluation. Compared 
to our previous definition, this definition significantly reduces the number of color bases 
for complicated processes. For instance, it is reduced from 24 to 14 1 for processes having 
two quark pairs and two gluons (e.g., gg — > qqqq). This reduction results in a faster 
calculation of the color factor and the color connection. 



iV jets (af ) 


base-subproc. 


subproc. w/ all flavors 


2 


8 


176 


3 


9 


276 


4 


14 


891 



Table 1: Number of base-subprocesses for the N jet production process in pp(p) collisions 
to be compared with the total number of subprocesses that they cover. The subprocesses 
can be classified according to the difference in the initial-state parton combination, and 
further classified by accounting for the charge and parity symmetries in the jet flavors. 



3 Program running 
3.1 Distribution package 

The distribution package of the GR@PPA event generator consists of two part, the frame- 
work and the matrix elements. The framework contains a CHANEL library, a BASES/ 
SPRING package, and a kinematics library. The CHANEL library is a set of subroutines 
to calculate diagram elements: vertices, propagators and external legs. BASES/SPRING 
is a multi-dimensional general-purpose Monte Carlo integration and event-generation pro- 
gram set. The kinematics library converts the random numbers given by BASES /SPRING 

1 Among 14 combinations, two color factors produce zero amplitudes at tree level. 
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to a set of kinematical variables for the matrix element calculation, and converts the re- 
turned matrix element to the differential cross section. These routines are common to 
all processes. On the other hand, the matrix elements, generated by GRACE with some 
modifications for hadron collisions, are process dependent. They are supplied as separate 
packages. 

The program running of GR@PPA generator is totally controlled by the subroutine 
GRCPYGEN. In the initialization, GRCPYGEN calls BASES to evaluate the total cross section 
for the given process, while it calls SPRING in the event generation cycle. The calling 
sequence of GRCPYGEN is as follows: 

CALL GRCPYGEN (CBEAM, IGSUB, MODE, SIGMA), 

where the input arguments are 



CBEAM (CHARACTER) 
IGSUB (INTEGER) 
MODE (INTEGER) 



'PP' for pp collisions and 'PAP' for pp collisions 
process number 

= 1 for calling BASES, and for calling SPRING, 



and the output is 

SIGMA (REAL*8) : integrated cross section. 

The argument CBEAM is dummy when MODE = 0. A unique number IGSUB is assigned to 
every physics process as listed in Appendix A. The output SIGMA is always equal to the 
integrated cross section of the process specified by IGSUB. The generated event is stored 
in the LHA common block jjj. 

The distribution package is arranged for the use on Unix systems. However, since the 
structure is rather simple, we expect that the program can be compiled and executed on 
other platforms without serious difficulties. The package for the GR@PPA framework is 
composed of the following files and directories: 



Conf ig.perl 

Makefile 

README 

VERSION-2.76 

proc . list 

basesv5 . 1/ 

chanel/ 

grckinem/ 

example/ 

inc/ 

lib/ 



a script to configure the setup, 

Makefile for the setup, 

a file describing how to set up the programs, 

a note for this version, 

a list of the processes supported in this version, 
BASES/SPRING (version 5.1) source codes, 
CHANEL source codes, 
source codes of kinematics, 
source codes of example programs, 
INCLUDE files, 

the directory to store object libraries; initially empty. 
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The proc .list contains all process names and the corresponding process identification 
numbers IGSUB. It also includes keywords for the processes used to configure the matrix 
elements. 

The matrix element packages are separately given. They need to be installed in the 
GR@PPA framework directory. We provide several packages. Users can choose any of 
them according to their needs. Each matrix element package is composed of the following 
directories: 

xxx : a set of matrix elements, 

diagram : Feynman diagrams used in the calculation of this process, 

where xxx is a process group name. For example, the matrix elements for top quark 
production processes (with or without extra jets) are all placed under the top directory. 

3.2 How to install 

The programs can be downloaded from 

URL: http://atlas.kek.jp/physics/nlo-wg/grappa.html . 



The first task for the setup is to run the configuration script Config.perl. Next, 
users have to edit the file Makefile to specify the paths to the libraries to be linked, 
such as PYTHIA, HERWIG and CERNLIB. Those parts to be edited can be found at the 
top of the Makefile. Users may also edit the inc/define.h file to select a PDF library. 
Currently, PDFLIB in CERNLIB |TS], LHAPDF PU, and PYTHIA built-in PDFs are 
supported. The standalone CTEQ6 [T7j routine is also included for the default use. In 
the case for using PDFLIB, edit one of the lines in inc/define.h as 



#def ine PDFLIB CERNPDFLIB . 



As an example we show an instruction to install the PU+ljet and Vy+2jets production 
processes in the GR@PPA framework. After downloading the framework package and the 
two matrix element packages, unpack them as 

tar xzvf GROPPA-2 . 76 . tgz 
cd GR@PPA-2.76 

tar xzvf . . /matrix_wlj_vl . 03 .tgz 
tar xzvf . . /matrix_w2j_vl . 03 . tgz . 

By these commands, the wjets directory is created and under it two directories, wlj and 
w2 j , are created. Then, configure the matrix elements and the compiling environment in 
Config.perl, and execute it as 

Config.perl . 



7 



This creates the Makefile according to the present matrix-element installation. Users 
have to edit the Makefile to specify further details, such as particular compiler options 
and paths to external libraries. Now, users can create the matrix element libraries and 
the framework libraries. Type the commands as follows: 

make wlj 
make w2j 
make kinem 
make integ . 

Note that, because the kinematics codes are modified according to the selected matrix 
elements, the kinematics library must be created after making the matrix element libraries. 
Users must repeat the above sequence again when they add new matrix elements. The 
created libraries are installed in the lib/ directory by executing the command, 

make install . 

The example programs are given for the standalone use, and PYTHIA- and HERWIG- 
interfacing. The command 

make example 

sets up these examples in the example/ directory. 

3.3 Initialization and customization 

Although the execution of GR@PPA is controlled by the subroutine GRCPYGEN, the de- 
tailed behavior depends on some parameters in common blocks and conditions defined 
in some subprograms. The subroutine GRCINIT initializes all related parameters to the 
default values. Users can change them after calling GRCINIT as described in the following. 

The parameter that is necessary to be given by users is GRCECM, which specifies the cm 
energy of the beam collision in GeV. Optionally, users can define some phase-space cuts 
in the laboratory frame: GPTCUT, GETACUT and GRCONCUT. These parameters define the 
minimum pt in GeV, the largest pseudorapidity in the absolute value and the minimum 
separation in AR, respectively. These cuts are applied to all produced jets except for 
those from weak boson decays. The separation (AR) is defined for every pair of jets as 

AR = ^Acj) 2 + At? 2 , (4) 

where A0 and A77 are the separation in the azimuthal angle and the pseudorapidity, 
respectively. The arrays GRCPTCUT, GRCETACUT and GRCRCONCUT can separately set certain 
phase-space cuts for each final state particle. 

Additionally, the subroutine GRCUSRCUT provides a framework to apply customized 
cuts referring to detailed properties such as four-momenta of the particles. Some useful 
functions are available there. If the events are acceptable IUSRCUT = should be returned, 
while IUSRCUT = 1 for the rejection. Since those cuts are applied at the integration stage, 
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a high event generation efficiency can be achieved in the event generation. The selection 
can be applied even for the particle flavor. This feature is useful for a parton flavor 
selection in the jet production. A concrete example can be found in Appendix B. 

The arrays IGWMOD and IGZMOD specify the decay mode of the W and Z bosons, 
respectively. If the value is equal to 1, the corresponding decay channel is activated; if 
the number is equal to 0, the channel is disabled. IWIDCOR is an option for the decay 
width correction for the W and Z bosons. The theoretically calculated values are used 
if IWIDCOR = 1, while the experimentally measured ones can be used if IWIDCOR = 2. If 
IWIDCOR = 2, the cross section is weighted by using the measured branching ratio as 

■ptot 

a(X -> //) = a calc -J* Br cxp (X -> //) , (5) 

calc 

where <r ca i c is the calculated cross section using the given total decay width r*°p, 
is the partial decay width theoretically calculated, and Br cxp (X — > //) is the branching 
ratio experimentally measured. Br exp (X — > //) must be given in the arrays GRCWBR and 
GRCZBR for the W and Z bosons, respectively. 

The array IGJFLV determines the active flavor of jets in the final state except for those 
from W and Z boson decays. IflGJFLV(I) = 1, the corresponding flavor I (d, u, s, c, 6, t, g 
for / = 1, 2, 7) is activated, while it is disabled if IGJFLV (I) = 0. Regardless that they 
are disabled by the IGJFLV flags, those flavors are occasionally produced if the diagram 
structure requires them. For example, consider the case that the u flavor in the PDF 
is activated and the qg — > qg base-subprocess is chosen in the QCD 2 jets production. 
In this case, if g is activated, the ug — > ug subprocess is activated even if u is disabled 
by IGJFLV. If both u and g are disabled, this subprocess itself is disabled. That is, the 
events are accepted if either of the activated flavors by IGJFLV has to appear in the final 
state. The jet flavor can also be controlled by changing the CKM parameter GRCCKM if 
the process includes intermediate W bosons. 

Further detailed parameters, such as particle masses, decay widths and couplings, are 
accessible in the subroutine SETMAS. The mass and the total decay width of the weak 
bosons can be manually controlled there. GR@PPA does not give any constraint to these 
parameters. However, the electroweak parameters are characterized by only three (plus 
Higgs mass) parameters. 2 The setting for dependent parameters are ignored and their 
values are recalculated according to the choice of the scheme given by the parameter 
IGAUGE. The scheme is taken as the default, where all the parameters are given by the 
set of following parameters: 

{G F , M w , M z ) . (6) 

Here, Gp is the Fermi constant, My/ and Mz are masses of W and Z bosons, respectively. 

The other process-specific parameters are defined in the subroutine GRCPAR. Users 
can choose different conditions for different processes. The variable ICOUP determines the 
energy scale (Q 2 ) for calculating the coupling strengths, a cm and a s , in the matrix element 

2 Additionally, a gauge invariance also demands fermion masses in the Yukawa sector. 
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calculation (renormalization scale). The selectable choices are prepared (see TableEJ). The 
variable IFACT determines Q 2 for PDF (factorization scale). The definition is the same as 
ICOUP. The same choice as ICOUP is taken if IFACT is not explicitly given. As an option, 
users can apply their own definitions of these energy scales, by setting ICOUP = 6 and/or 
IFACT = 6 and editing the subroutine GRCUSRSETQ. An example is attached to grcpar.F. 
The number of flavors used in the coupling calculation and PDF is given in INPFL. The 
flavor is ordered as u, d, c, s and b. The extrapolation to the top quark in the coupling 
formula and PDF is not taken into account at the moment. 

The integer variable IBSWRT controls whether BASES should be called in the initial- 
ization or not. The task of BASES is to optimize the integration grids and, after that, 
store the optimized results in a "BASES table" . The execution of BASES consumes much 
CPU time because a precise evaluation is necessary for an efficient event generation by 
SPRING. It is not necessary to repeat the execution for identical conditions. A previously 
optimized result ("BASES table") is reused if IBSWRT = 1. It should be noted that, once 
the parameters which affects the estimated result of the integration are changed, the "con- 
dition" is no longer identical and BASES has to be re-executed. The parameter NCALL 
specifies the number of sampling points in each step of the iterative grid optimization in 
BASES. The larger this number is, the better the conversion would be. However, it takes 
longer in the CPU time. The optimized values are preset in grcpar.F. The character 
variable GRCFILE gives the "BASES table" file name 3 . A new file must be specified if 
IBSWRT = 0, while an existing file must be specified if IBSWRT = 1. 

The parameters described in this subsection are summarized in Table El 

3.4 BASES integration 

In the output of GR@PPA, users should pay appropriate attention to the print out from 
BASES, especially when they apply tight cuts. Since each subprocess is composed of 
many coherent diagrams, it is not practicable to take all singularities into account in the 
"kinematics" definitions. Some very minor ones are ignored in GR@PPA. A combination 
of very tight cuts may enhance the relative contribution of ignored singularities. In such 
cases, it is likely to happen that, in the BASES iteration, the estimated total cross section 
jumps (increases) to a value unreasonably different from the previous estimation and, 
accordingly, the estimated error increases. Users should consider that they must be in 
such a trouble if they find a jump of, for instance, more than three times the previous 
error. The results are unreliable in the phase-space region defined by such cuts. The 
instructive integration accuracy is 0.5% or better for every iteration. Users should change 
the parameter NCALL to a larger value if this accuracy is not achieved. 

3 BASES actually creates two files having extensions of .data and .result, respectively, added to the 
name given by GRCFILE. The former is the "BASES table" , while the latter is a readable summary of the 
BASES execution. 
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Parameter 


Description 


GRCECM 


CM energy of the beam collision in GeV. (D=14000.) 


CBEAM 


'PP' for pp collisions and 'PAP' for pp collisions. (D='PP') 


IGSUB 


Process number. See Table 13 


MODE 


Mode selection for GRCPYGEN: = 1 for calling BASES, and for 
calling orrvllNLj. 


GPTCUT 
GRCPTCUTC8 - ) 


Minimum pt cut in GeV tor jets except tor those from weak 
boson decays (D=20 OeVl 

TYTinimnm Tim put, in CtpV for p^ph fin^l ^tatp n^rtiplp |T^=D CtpA/i 

IV llllllllLllll LJ 1 V LX U 11X V — 1 \j V X V./X v_<CX v^XX XXXXCXX O U CX UVj IJCXiX U X V . XV . 1 J / \J . V_. 1 V_, V / 


GETACUT 
GRCETACUTC8) 


1 j^TP'P^it, ■nspnrlnTRnirli'hv put in thp Rhsnhit.p vrIup tor ipt^i p^ppnt, tor 

1 J CXiX CVjO U UOvj LiVXvl CXi l_XXV_X.X I V v. LX U XXX U llVj CXi UOV/1 \JL UVj V CXjX LX vj 1U1 1 Vj U O vAvb l_X L) X V_y X 

those from weak boson decays. (D=3.) 

Farthest nsendoraniditv cut for each final state narticle (D=10 ) 

l J CAjX V_jkJ> U IJ kjl LivXv/1 LA l_y 1V11 U y Vj I. i. u Ivl LX X X XXXXtXJX kjj U CXi U l_y CXX U X v . L * \ J s X- VX * J 




Minimum separation cut in Ai? (= v/ A</> 2 + A?] 2 ) for jets except 
for those from weak boson decays. (D=0. 4) 


GRCRC0NCUT(8) 


Minimum separation cut in AR (= a/A0 2 + A77 2 ) for each final 
state particle. (D=0.) 


IGWM0DC20) 


Decay mode for W boson: = 1 to activate, and to deactivate. 


IGZM0DC16) 


Decay mode for Z boson : = 1 to activate, and to deactivate. 


IWIDCOR 


Width correction for W and Z. (D=l) 


GRCWBRC20) 


Branching ratio tor W . ft is used only when IWIDC0R=2. 


GRCZBRC16) 


Branching ratio for Z. ft is used only when IWIDC0R=2. 


IGJFLVC7) 


Flag for the jet flavor of jets except for those from weak boson 
decays: = 1 to activate, and to deactivate. 


IGRCGEF 


Photon interference in Z production processes is included 
if set to 1, ignored if 0. (D=l) 


GRCCKM(3,3) 


CKM parameters. 


IGAUGE 


Choice of electroweak parameters. 

The scheme is used as the default. (D=l) 


I COUP 


Choice of the renormalization scale. 

= 1 : \f§ of the hard interaction. 

= 2 : average of squared transverse mass (< m\ >). 

= 3 : total squared transverse mass (Yl m\). 

= 4 : maximum squared transverse mass (max m\). 

= 5 : fixed value. Set GRCQ in GeV. 

= 6 : user defined scale. Set GRCQ in the subroutine GRCUSRSETQ. 


I FACT 


Choice of the factorization scale. (D=0) 
The definition is the same as ICOUP. 

If IFACT = 0, the same value as the renormalization scale is used. 
In case IFACT = 5 or 6, set GRCFAQ in GeV. 


GRCFILE 


Output file name for the BASES integration. 


IBSWRT 


Mode selection for BASES integration: = for calling the BASES 

integration, and 1 for skipping. (D=0) 

If IBSWRT = 1, the file defined in GRCFILE is used. 


NCALL 


IN umber of sampling points in each step of the iterative grid 
optimization in BASES. 


INPFL 


Number of flavors used in the coupling calculation and PDF. (D=5) 



Table 2: Parameters in GR@PPA. 
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4 Processes 



4.1 V + N jets 

The process identification number IGSUB is assigned to be 100 ~ 104 and 110 ~ 114 for 
W and Z production processes with up to 4 jets. The matrix elements (ME) include the 
boson decays into two fermions, so that the decay properties are correctly reproduced at 
the tree level. Heavy quarks such as bottom and top quarks can also be included as the jets 
as well as the lighter quarks and the gluon. Since the quark mixing in the couplings with 
W boson is treated up to the third generation, single and triple heavy quark productions 
are allowed in the multi-jet configuration, in addition to the pair production from the 
gluon splitting. In Z + N jets processes, the Z/Y interference can be taken into account 
in the ME calculations. The interference between decay fermions and those from other 
sources is ignored. 

The benchmark cross sections of the V + N jets production processes for the Tevatron 
Run II and the LHC conditions are presented in Table EJ where the quark flavors are 
included up to the b quark. We also present the benchmark cross sections for the case 
having at least one b quark, and having a tt pair in the final state. The results are in good 
agreement with those from other generators [TS1 IP*?! HU\ |2~T] 12*2] . The detailed parameters 
used in these calculations are described in Appendix B. 





Tevatron Kun-11 


LHC 


iV jets 


W(ev e ) 


Z(e+e-) 


W(ev e ) 


Z(e+e-) 



1 
2 
3 
4 


1.576(2) xlO 3 
1.852(3) xlO 2 
3.461(7) xlO 1 
6.29(2) 
1.201(5) 


1.598(3) x 10^ 
1.829(4) xlO 1 

3.485(6) 
6.35(2)xl0- 1 
1.173(3)xl0" 1 


1.116(2)xl0 4 
2.854(5) xlO 3 
1.143(3) xlO 3 
4.82(l)xl0 2 
2.19(l)xl0 2 


9.57(3) xlO 2 
2.614(7) xlO 2 
1.082(2) xlO 2 
4.53(l)xl0 1 
2.045(5) xlO 1 


2 (> 1 b) 

3 (> 1 b) 

4 (> 1 b) 


3.260(6)xl0" 1 
1.019(2)xlO _:L 
2.947(8) xlO" 2 


9.24(2) x 10^ 
2.266(5)xl0~ 2 
3.817(5)xl0" 3 


2.720(7) 
4.305(9) 
3.90(3) 


8.68(2) 
3.740(7) 
9.47(2) xlO^ 1 


tt + 
tt + 1 
tt + 2 


5.269(9) xlO" 4 
1.357(2)xl0" 4 
6.86(l)xl0- 5 


2.248(3) xlO" 4 
6.302(9)xl0- 5 
5.707(5) xlO" 7 


3.774(7) xlO" 2 
4.98(2) xlO" 2 
6.52(2) xlO" 2 


2.682(6) xlO" 2 
3.59(l)xl0- 2 
1.156(2)xl0" 2 



Table 3: Benchmark cross section (pb) for V + iV jets processes. Results are presented for 
the Tevatron Run II and the LHC cases. The detailed parameters used in the calculations 
are described in Appendix B. 

4.2 VV + N jets 

The number IGSUB is assigned to be 120 ~ 122 for the WW production processes, 130 ~ 
132 for the WZ production processes, and 140 ~ 142 for the ZZ production processes. 
The matrix elements (ME) include the boson decays into two fermions, so that the decay 
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properties such as the spin correlation between two bosons is correctly reproduced at the 
tree level, while the interference between the fermions from different bosons and those 
from other sources is ignored. The Z /j* interference can be taken into account in the 
ME calculations. Heavy quarks (b and t) can also be included in the jets. 

The benchmark cross sections for the Tevatron Run II and the LHC conditions are 
presented in Table 0] We also present the results for the case having at least one b quark, 
and having a tt pair in the final state. The detailed parameters are described in Appendix 
B. 





Tevatron Kun-11 


7V jets 


WW 


WZ 


ZZ 



1 
2 


7.91(2)xHT* 
1.986(5) xl(T 2 
4.77(l)xl0~ 3 


6.06(l)xl0~ 3 
2.013(4) xl0~ 3 
5.16(2)xl0" 4 


1.541(3) xlO" 3 
3.721(6) xlO" 4 
7.698(9) xlO" 5 


2 (> 1 b) 


1.582(2) xlO" 4 


6.00(2)xl0- b 


1.629(3) xl0- b 


tt + 


5.504(5) xlO" Y 


5.065(7) xl(T 8 


8.871(9)xl0- y 






LHC 


7V jets 


WW 


WZ 


ZZ 



1 
2 


1.339(5) 
1.071(3) 
9.41(2) xlO -1 


4.64(l)xl0~ 2 
1.873(4) xlO- 1 
1.536(4) xlO" 1 


1.080(2) xlO -2 
6.06(l)xl0~ 3 
3.419(5)xl0" 3 


2 (> 1 b) 


7.36(l)xl0" 2 


1.164(3)xl(T 4 


1.887(8) xlO" 4 


tt + 


7.432(9) xlO" 4 


1.213(2)xl0~ 5 


3.040(5) xl0- e 



Table 4: Benchmark cross section (pb) for VV + iV jets processes. Results are pre- 
sented for the Tevatron Run II and the LHC cases. The detailed parameters used in the 
calculations are described in Appendix B. 

4.3 bbbb 

We provide the bbbb production processes which have been developed in our previous work 
[TUj as a separate matrix element package. The process identification number is 160 for 
the Higgs production associated with bottom quarks (Y" b 2 a 2 ), 161 for the Z/j* mediated 
process (a 2 ttg m ), 162 for the pure QCD process (a 4 ), 163 for the HZ production (F 6 2 a 2 m ), 
and 164 for the Z pair mediated process (a 4 m ), where the SM Higgs boson is assumed for 
the Higgs production. 

4.4 tt + N jets 

Aside from the QCD multi-jets processes, we provide top pair production processes. In 
these processes, the whole decay chain of the top quark is included in the diagram cal- 
culation, so that the spin correlation in the top decay is fully reproduced. In addition 
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to the tt production (6-body ME), we also provide the tt + 1 jet (7-body ME) process. 
This process may give an insight for the 7-body kinematics including the top decay. The 
process number is assigned to be 170 for tt, and 171 for tt + 1 jet. Possible anomalous 
couplings in the top decay and its production will be included in the ME calculation in a 
future version. They will be provided as an update set of the matrix elements. 

4.5 N jets 

The QCD multi-jets processes are provided with the process identification numbers from 
182 to 184 for 2, 3 and 4 jets, respectively. The electroweak interaction is ignored in these 
processes. Heavy flavors such as bottom and top quarks can be produced by appropriately 
setting IGJFLV. 

The benchmark cross sections of the QCD N jets processes are presented in Table 03 
We also present the results for heavy flavor production processes. These results are in 
good agreement with those from other generators [THl EH EU 122] • It should be noted 
that the result for the bbbb production is the same as that for the process IGSUB = 162. 
The detailed parameter setting is described in Appendix B. 





Tevatron Kun-11 


LHC 


N jets 


QCD jets «) 


QCD jets (af ) 


2 


1.853(5)xlO Y 


4.355(8)xl0 8 


3 


6.88(3)xl0 5 


3.326(7) xlO 7 


4 


8.82(3)xl0 4 


7.94(1) xlO 6 


bbbb 


5.70(1) 


8.42(2) x 10* 


bbtt 


7.46(2) xl0~ 3 


3.796(8) 


tttt 


5.809(4) xlO" 6 


2.322(5) xlO" 2 



Table 5: Benchmark cross section (pb) for QCD multi-jets processes. Results are pre- 
sented for the Tevatron Run-II and the LHC cases. The detailed parameters used in the 
calculations are described in Appendix B. 



5 Performance 

The computation performance of GR@PPA for the W + N jets processes in the Tevatron 
Run-II condition is summarized in Table |3 The tests have been performed using Intel 
Pentium 4 3.4 GHz CPU and two different Fortran compilers: a free software, g77 in gcc 
2.96, and a commercial compiler, Intel Fortran Compiler version 8.0. The integration time 
and the generation speed are separately shown. Clearly, the commercial compiler is about 
2.5 times faster than the free compiler. However, in both cases, the time is not intolerable 
for large scale Monte Carlo productions. The generation efficiencies of SPRING, shown 
in the table, are exceptionally good for this kind of complicated processes. 



14 



The result for the W + 4 jets is not shown in the table, since it consumes too long 
CPU time to run as a single process. A parallel computing is indispensable for practical 
uses. A parallel computing employing MPI (Message Passing Interface) is supported 
in BASES. The benchmark result in Table El has been obtained by using a PC farm 
having 12 Intel Pentium 4 3.2 GHz CPUs with Intel Fortran Compiler version 7.0. The 
integration time is 23 hours fully using these CPUs. The event generation speed is 0.11 
events/sec with the generation efficiency of 0.1%. Though SPRING does not support 
parallel computing, users can parallelize the event generation by hand by changing the 
starting random number seed. The parallel computing has also been applied to the W + 
3 jets production. The integration time has been reduced to 24 minutes from 5 hours for 
the single-process computing on the same platform. The scalability is quite good. Those 
users who want to apply the parallel computing should contact the author 4 . 



Process 


Fortran 


Int 


egration 


Event Generation 


Efficiency 




Compiler 


time 


(H:M:Sec) 


speed (events/sec) 


(%) 


W(ev e ) + jet 


g77 


00 


00 


05 


43859 


69.6 


intel 8.0 


00 


00 


02 


101010 




W{ev e ) + 1 jet 


g77 


00 


00 


52 


13927 


19.9 


intel 8.0 


00 


00 


19 


34364 




W{eu e ) + 2 jets 


g77 


00 


38 


27 


709 


1.7 


intel 8.0 


00 


13 


52 


1957 




W(eu e ) + 3 jets 


g77 


14 


03 


46 


18 


0.3 


intel 8.0 


04 


57 


50 


52 





Table 6: Performance of GR@PPA for the W + N jets processes in the Tevatron Run-II 
condition. The tests have been performed using Intel Pentium 4 3.4 GHz CPU and two 
different Fortran compilers: g77 version 2.96, and Intel Fortran Compiler version 8.0. The 
integration time and the generation speed are separately shown. 

6 Summary 

The GR@PPA event generator has been updated to version 2.7. In this version, the 
extensions have been mostly applied to deal with the flavor configuration in the initial 
and final state by sharing several subprocesses into the single subprocess. This allows us 
to incorporate the variation in the initial and final states parton configurations such as 
the multi-jets production processes. A new method has been also adopted to calculate 
the color factors of complicated QCD processes. Those implementations are able to speed 
up the diagram calculation for multi-jet production processes significantly. 

The distribution of this version is composed of two separate packages: the frame- 
work and the matrix element for the processes. The framework consists of some process- 
independent elements such as integration packages and kinematic libraries. The matrix 

4 Soushi.Tsuno@cern.ch 
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element provides a set of the matrix elements used in the GR@PPA event generator. 
Currently, we provide the matrix element packages for V (W or Z) + jets (< 4 jets), 
VV + jets (< 2 jets) and QCD multi-jet (< 4 jets) production processes at the tree level 
for pp and pp collisions, where the jet flavor includes up to top quark. The four bottom 
quark productions implemented in our previous work (GR@PPA_4b) are also included. 
In addition, we provide the top-pair and top-pair + jet production processes, where the 
correlation between the decay products are fully reproduced at the tree level. Namely, 
processes up to seven-body productions can be simulated, based on ordinary Feynman 
diagram calculations. 

The GR@PPA event generator takes advantages that further extensions are easily 
applicable in the GR@PPA framework. Once NLO processes or non-Standard Model 
processes are given as a matrix element package, users will be able to simply use them 
without any detailed care. An extended framework to include the parton shower is also 
possible. They may be available in a future release or by a release of the new matrix 
element packages. The GR@PPA event generator will be suitable not only for a large 
scale Monte Carlo production for high luminosity hadron collisions at Tevatron and LHC, 
but also for future NLO calculations to be composed of lots of subprocesses. 
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A Processes in GR@PPA 

All physics processes are recognized by unique process numbers IGSUB listed in Table 
The second column in the table is the order of coupling constants for the process. The 
parameter is the Yukawa coupling of the b quark. The third column shows the command 
paramter for compiling the ME. The parameter name is the same as the directory name 
of the process. The forth column is the process descriptions. Each process consists of 
several base-subprocesses. The symbol "/" denotes the fermion from boson decays. The 
decay fermions except for those in the bbbb processes do not interfere with other partons. 
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IGSUB 


Coupling 
order 


Command 
parameter 


Process description 


100 
101 
102 
103 
104 


a 1 

2 

a em a s 

a 2 a 2 

2 3 

or or 

9 A 

a a 


wOj 
wlj 
w2j 
w3j 
w4j 


pp(p) - W(2f) + jet (+X) 
pp(p) -> W(2/) + 1 jet (+X) 
pp(p) - W(2/) + 2 jets (+X) 
pp(p) - W(2/) + 3 jets (+X) 
ppip) -> W(2/) + 4 jets (+X) 


110 
111 
112 
113 

114 


a 2 

^ cm 

a 2 a 2 

em s 

a 2 a 

^2 ^4 
ft ft 
"em u s 


zOj 

zlj 
z2j 

z3j 

z4j 


pp(p)-Z/ 7 *(2/)+0jet (+X) 
pp(p) - Z/ 7 *(2/) + 1 jet (+X) 
pp(p) -Z/ 7 *(2/) + 2 jets (+X) 
pp(p) -Z/ 7 *(2/) + 3 jets (+X) 
-> Z/7 (2/) + 4 jets (+A) 


120 
Izl 
122 


ft 4 

"cm 

4 

ft 4 ft 2 


wwOj 
wwl j 
ww2j 


pp(p) - W^(2/)^-(2f) + jet (+X) 
PPIP) H / ^(2/jl4/ (2/ j + 1 jet (+AJ 
pp(p) -> W" + (2/)iy-(2/') + 2 jets (+X) 


130 
lol 
132 


ft 4 

cm 

4 

ft cm fts 

ft 4 ft 2 
em s 


zwOj 
zwlj 
zw2j 


pp(p) - Z/Y(2f)W(2f) + jet (+X) 
PPVP) Z/ 7 (2/ J 14/(2/ J + 1 jet (+AJ 
pp(p) - Z/ 7 *(2/)W(2/') + 2 jets (+X) 


140 
141 
142 


ft 4 

"cm 
4 

ft ft 
em 5 

ft 4 ft 2 


zzOj 
zzlj 
zz2j 


ppip) - Z/ 7 *(2/)Z/ 7 *(2/') + jet (+X) 
pp{p) - Z/ 7 *(2/)Z/ 7 *(2/0 + 1 jet (+X) 
- Zh*(2f)Zh*(2f) + 2 jets (+A~) 


160 
161 
162 
163 
164 


a 2 a 2 

ft 4 

o "em 

ft 4 

cm 


b4hbb 
b4zbb 
b4qcd 
b4hz 
b4zz 


-> h (bb) + 66 (+A~) 
pp(p) -> Z/ 7 * (66) + 66 (+X) 
pp(p) -> 6666 (+X) 
#p(p) -> ^o(66) + Z/ 7 * (66) (+X) 
pp(p) -> Z/ 7 *(66) + Z/ 7 * (66) (+X) 


170 
171 


ft 4 ft 2 
em s 

Ot or 


tt6bdy 
ttj7bdy 


pp(p) -> tt -> 6/ (+X) 
pp(p) -> tt + 1 jet -> 6/ + 1 jet (+X) 


182 
183 
184 


ft 1 


qcd2j 
qcd3j 
qcd4j 


pp(p) — > 2 jets (+A~) 
~ 3 jets (+A~) 
ppip) -> 4 jets (+X) 



Table 7: Processes included in GR@PPA, version 2.7. 
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B Parameters 

The kinematical cuts and other parameters used in the benchmark tests in Section 4 are 
as follows: 

Beams : 

The benchmark tests have been carried out for the Tevatron Run II (pp) and the LHC 
(pp) conditions with the center-of-mass energies of 

1.96 TeV pp (Tevatron Run II), , , 

14 TeV pp (LHC). ^ ' 

Electroweak parameters : 

The scheme is adopted with the parameter set of 

(G F , M w , M z ) = (1.16639 x 1(T 5 GeV" 2 , 80.419 GeV, 91.188 GeV) , (8) 



where Gp is the Fermi constant, Mw [23] and Mz [23] are the masses of the W and Z 
bosons, respectively. This set leads to the other parameters as 

a em = 1/132.51 , sin 2 ^ = 0.2222 . 

Boson width and decay : 

The W boson is forced to decay to a pair of an electron and a neutrino, and the Z 
boson to an electron-positron pair. In diboson production processes, the interference 
between electrons (neutrinos) from two bosons is ignored. We use fixed decay widths 
for W, Z and the top quark as 



w 



2.048 GeV, T z = 2.446 GeV, T t = 1.508 GeV. (9) 



CKM parameters : 

We use the following CKM parameters: 



d' \ 


. ( \v ud 


| = 0.9752 | 


v us 


| = 0.2210 | 


v ub 


| = 0.0054 


s' 


= WJ 


= 0.2210 


v cs \ 


= 0.9743 


\v cb \ 


= 0.0419 


«) 


\ \Vu\ 


= 0.0054 


\Vu\ 


= 0.0419 


\Vt b \ 


= 0.9991 




(10) 



The phase of CKM matrix is neglected. 

Fermion masses : 

Following values are use for the masses of heavy quarks 

m c = 1.5 GeV m b = 4.7 GeV m t = 174.3 GeV . (11) 

The other quarks (u,d,s) are assumed to be massless. The electron and neutrino masses 
are also neglected. 
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Proton Distribution Function : 

We use CTEQ6L ^7] containing 5 flavors for PDF. The corresponding strong coupling 
constant is a s {M\) = 0.1180. 

Energy scale : 

The renormalization and factorization scales (Q) are always chosen to be identical and 
fixed to the Z boson mass as 

Q — IMf — fJ,R — mz — 91.188 GeV. (12) 

Kinematical cuts : 

The following kinematical cuts are applied to avoid infrared/collinear singularities, 

p T > 20 GeV , \r]\ < 3.0 , AR > 0.4 . (13) 

These cuts are applied to the jets (u, d, c, s, b and g) and electrons, while no cut is 
done to top quarks and neutrinos. 

Flavor selection : 

In some tests of heavy flavor productions, a particular jet flavor is selected in the 
integration. This is done in the subroutine GRCUSRCUT. The following is an example to 
select those events having two top quarks in the final state: 

IT0P = ! counter of the number of top-quark in jets. 

DO I = l.NJET ! loop over n jets. 

IF (IABS (GRCKFCD (I) ) . EQ . 6) IT0P = IT0P + 1 ! number of top qaurks . 
ENDD0 

IF (IT0P.NE.2) IUSRCUT = 1 ! reject this event. 
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